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Abstract 

Using a Hagedorn resonance gas picture and quark-hadron duality we estimate the dilepton 
emission rate in the vicinity of the QCD deconfinement phase transition. The result is then 
used to calculate a dilepton spectrum in ultrarelativistic heavy ion collisions. We show that 
multibody contributions taken into account in the Hagedorn resonance gas approach provide 
an enhancement of the production rate of massive dileptons as compared to the previously 
considered sources. 
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1 Thermal rate from a dense hadron matter 



The dilepton production in high energy heavy ion colhsions has drawn considerable attention 
as a possible signal for the appearance of a quark-gluon plasma see also the reviews [§. 
The main problem, however, is to distinguish the signals coming from the normal hadronic 
phase and from the quark-gluon plasma. To be able to do so one needs reliable models for the 
dilepton production. At extremely high temperatures the situation is relatively clear. Due 
to the asymptotic freedom of QCD the dominant production mechanisms are perturbative. 
At very low temperatures we also have a reliable description of dilepton production in a 
rarefied hadron gas, predominantly by pion annihilation. Much less can be said about the 
intermediate temperature range around the phase transition temperature Tc ~ 160 MeV, 
where the hadronic gas is very hot and dense and the usual methods do not work, but the 
quark-gluon plasma is cold and rarefied and the perturbative QCD results are less reliable. 
We shall address the problem of estimating the dilepton emission from a hadron gas at these 
intermediate temperatures. This is of special interest because a hadron gas close to the 
transition temperature Tc could be a dominant feature of the final state produced at the 
present CERN heavy ion experiments with lead beams. 

The new experimental data on dilepton production in heavy ion collisions demonstrate 
a pronounced excess of produced dileptons over those produced by known hadron sources 
both in low-mass range and at larger masses 0]. This is a clear motivation to look for 
additional sources of dileptons not taken into account in the conventional treatment. The 
excess of small- mass dileptons has been extensively discussed (see, e.g., H). More related 
to our discussion is, however, the excess in the intermediate mass range 1.5 GeV < M < 2.5 
GeV by roughly a factor of 2 with respect to expected contributions from the known sources. 
This is because the theoretical description we are developing is valid only for sufficiently large 
invariant masses of lepton pairs. The results will be shown for the mass interval 1.5 GeV 
< M < 4 GeV. We shall argue, that by improving the treatment of dilepton emission from 
the dense and hot hadron gas at the vicinity of one can get such an increase with respect 
to the conventional description. 

In we have proposed to use a Hagedorn resonance gas picture and quark-hadron 
duality to calculate the dilepton emission rate in the hadron phase at temperatures close to 
the critical one. Below we shall refine these arguments and present additional support to the 
main conclusion of . In the next section we shall construct an interpolation for the dilepton 
emission rate satisfying the conditions following from the quark - hadron duality and use 
it in a hydrodynamical model for heavy ion collisions. In the last section we formulate our 
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conclusions. 

The problem of calculating the dilepton emission rate in a hadron gas including also the 
mass range 1.5 GeV < M < 4 GeV has been considered in H and M. We compare the 



results of |q, 10 1 for the rate from binary hadron collisions with our results from the dense 



resonance gas at T^. We shall comment on the relation of our results to those of 0. Since our 
treatment effectively averages over mass intervals of the order of vector meson widths and 
mass differences, we cannot address the problem of medium effects on individual resonance 
parameters, relevant for the low-mass region of CERES results 0]. 

Let us first remind of the picture of dilepton production at a critical temperature proposed 
in 0. The main idea can be formulated as follows. In the vicinity of the critical point 
the hadronic matter is extremely dense B and the interactions between the hadrons are 
very strong. The idea of a Hagedorn resonance gas approach [|T^] is to encode the infinitely 
complicated strong interaction physics of a hot and dense hadron matter in a single function, 
an exponentially rising mass spectrum 

p(m) = c— e'"/'^° . (1) 

The main idea proposed in [0 was that the exponential spectrum in Eq. (1) takes care 
of strong interactions only. The free massive degrees of freedom saturating the partition 
function of the Hagedorn gas can still experience electromagnetic and weak decays. Of a 
particular interest to us is the spectrum of vector mesons, which are the only important 
source of dilepton production in the Hagedorn resonance gas in the considered mass range. 

In the vector dominance model, to be used below, the dilepton width of a vector meson 
is given by 

_ ^_4a^ ... 
~ g\M) 3 ■ ^ ' 

The spectrum of dileptons produced by the decays of vector resonances belonging to the 

Hagedorn gas can be written in the relativistic kinetic theory in the form 

= J ^"^^^("^^ J 2E(2^ -M-E/T) J 2^^(2.)3 2i,,(2vr)3 
X \Mip^pi+p2W (27r)45W (p-Pi- P2mpi + P2f - , (3) 



where p is the four-momentum of a decaying vector meson with a mass m, E = -^/p^ + 
is its energy, 2 are the four-momenta of the emitted lepton and antilepton, and py(M) is 
the subspectrum of vector mesons in Eq. (|^), not necesserily exponential. In Eq. (||) we have 
assumed a Boltzmann distribution for the energy of the decaying meson which holds with a 
very good accuracy provided the invariant mass is high enough. 
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In the above formula the squared matrix element for the V l^l decay, including spin 
summation and averaging, has the form 



Ana \ A , 



\g{m)) 3 

which also gives the the result Eq. (j^) for the decay width. For the dilepton emission rate 
one gets from Eq. 

The most important point is now to fix the mass dependence of the coupling constant 
1/ g{M). Here we can use the known cross section of e^e~ annihilation into hadrons which 
takes place in the same vector channel we are considering. The cross section of e^e~ anni- 
hilation into heavy vector mesons characterized by the mass spectrum py{M) reads: 

where Q, the CMS energy of e"'"e~ collision, equals the mass of produced vector meson. 
We proceed by assuming in the spirit of the Hagedorn bootstrap picture that the total 
cross section of e^e~ annihilation into hadrons is saturated by the production of vector 
meson states. This allows us to relate the emission rate from the Hagedorn gas to the total 
e+e^ hadrons cross section. We first have 

Ait cy^ 

c,(e+e- ^ hadrons) = i?-P(g) — — = a(e+e- V){Q) . (6) 

3 

From the above two equations we get the the mass dependence of the coupling constant 

1 _ i?^^P(M) 1 1 
g\M) ~ 27r2 M pv{M) 



(7) 



and finally, by inserting this into Eq. (|^), we can calculate the dilepton production rate in 
terms of known quantities: 



jpi+i- 2 



Let us rewrite the function R^^'^iM^ as 



R^^P(M) = R^"'-^ + bR{M) , (9) 
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where we have decomposed R'^^^{M) into a sum of its parton model prediction Rp°''^^ = 
NcJ2f ^"jy where e/ is an electric charge of the quark in units of the electron charge, and the 
deviation from the parton model result, 6R{M). Experimentally one has R^^^ ~ j^part ^j^]-^ 
a good accuracy for the mass-averaged cross section even though the local mass spectrum 
shows a visible resonance structure parametrized by 6R. From Eqs. (||,|) we get 

^f?i+i- ^2 2 

-JTTT = R'"'''^MTK^{M/T) + 6R{M)—MTK,{M/T), 

or finally 

This result shows that the parton model emission rate of dileptons with mass M equals that 
from the resonance gas at T ~ with the same accuracy as parton model describes the total 
cross section for e"'"e~ hadrons at the CMS energy M. We should like to stress, that this 
result does not depend on a particular assumption on the form of the heavy vector meson 
spectrum. A natural interpretation of the result is that hadrons as local annihilating degrees 
of freedom do not provide a simple description of the dilepton production in the immediate 
vicinity of a critical point, where they start to dissolve to their partonic constituents. As the 
correlations among the constituents get weaker the electromagnetic current, when probed 
through the production of heavy lepton pairs, is effectively that of quarks. 

This picture leads us in turn to a conjecture on the general behaviour of the dilepton 
rate as a function of temperature. Let us recall, that a formal nonperturbative expression 
for the dilepton rate in the hadron gas reads 

dN o? 



d^q Gvr^g^ 



J d'^xe-'"'' « J^{x)Jf,{0) »T . (11) 



Basing on Eq. (|T^) and the above- described physical picture it is tempting to suggest that 
at high enough dilepton masses the current- current correlator in (^) and thus the dilepton 
emission rate is to a good accuracy a continuous function of temperature across the phase 
transition point. How low in the dilepton invariant mass does this hold is determined by the 
accuracy related to the approximation made in Eq. (^. It would be interesting to check this 
conjecture in lattice stimulations. 

As our aim is to calculate dilepton production in heavy ion collisions, we shall need 
a continuous rate interpolation from the hadron rate to the QCD one as a function of 
temperature. A simple way of constructing an interpolation which will agree with Eq. (|10D 
is to start from the rate due to binary hadron collisions calculated in ^ and match it to the 
QCD one at the critical temperature T^. 
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Before turning to the actual hydrodynamic calculation, let us also mention a calculation 
in based on using the soft pion theorems in calculating the hadronic contribution to the 
polarization tensor. The main result of this paper for thermal dilepton rate from the hadron 
gas can be rewritten in the form 



dM^ dM^ 



;i + SR/W"'') , (12) 



2 p'^'^iM) 

where e = T^/6F^, = 93 MeV and p*^"^, p^ and p^ are the suitably defined electromag- 
netic, vector and axial spectral densities which can be extracted from measured quantities 
(for details see [^]). In the following we shall be interested in the production of dilepton pairs 
having invariant masses M > 1.5 GeV. In this mass range one can neglect p^ with respect 
to p^ and take p^ = p^"^ 0. This gives, when neglecting Ji?, 

^ - _ (, _ £!), (13) 

dhP ' dAP ' ^ 2"- ' ' 

It is interesting that the temperature dependence through e is in fact quite smooth, so the 
rate is predicted to be close to the thermal quark- ant iquark one even at temperatures below 



Tc. This obviously agrees with our main result Eq. (10). Unfortunately the range of masses 
where experimental information is available for the technique of is limited to M < 2.2 
GeV, so to analyze a broader mass interval one in forced to use a more phenomenological 
approach as described below. 

Let us now turn to the calculation of the dilepton emission rate from a hadron gas 
due to binary collisions of mesons ||^. In view of the previous discussion it is of special 
interest and importance to look how close is the rate of this approach to our result. This 
will provide valuable information on the relative importance of the higher resonances and 
the multibody effects in dilepton emission processes. In the Boltzmann approximation for 
the thermal distribution, which works fairly well when the invariant masses of lepton pairs 
satisfy M > 1.5 GeV, the rate reads 

'-^--'"i^mM/T). (14) 

where ah{M) is calculated from the specific Lagrangian for meson interactions [§. In the 
following we shall use the rate Eq. ( p!4|) including the following reactions Q (for brevity we list 
only the incoming particles) 0: tt+tt", K+R-, K^K^, p+p-, K*+K*~, K*^K*^, up, tx^uj, 
■n^cj), r]p^, r]uj, r](f), r/'p", r/'cj, r/'0, K+K*- + K-K*+, K^K*^ + K^K*^, 7r°u;(1420), 7r°u;(1660). 

^We are very grateful to P. Lichard for providing us the numerical data for this rate. 
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Figure 1: Ratio of static hadronic dilepton production rate in a binary collision approxima- 
tion to the quark- antiquark one at T=160 MeV 

In Fig. |I| we show the ratio of this binary hadron rate to our result Eqs. and ( ]T(]| ) at 
T = 160 MeV in the mass interval of interest. We see that the rate from binary interactions 
is less than our rate which takes into account also higher resonances and multibody effects 
and, in the approximation we use, equals the rate from quark- antiquark annihilations at this 
temperature. In the Hagedorn resonance gas approach the critical temperature is expected 
to be of order (Tc ~ 160 MeV), so taking into account contributions from binary collisions 
only is not sufficient at least in the vicinity of the critical point. In the end of the next 
section we shall discuss the numerical importance of this in calculating dilepton production 
in heavy ion collision in a hydrodynamic model, but already now one could notice, that the 
dependence on temperature in the above ratio cancels (cf. Eqs. and (0), so the difference 
is directly translated into the changes in the total dilepton rate. 

2 Dilepton emission in a hydrodynamic model 

In this section we shall study the implications of the above considerations on dilepton pro- 
duction in heavy ion collisions described with the simplest one-dimensional scaling hydrody- 
namic model [|15|. The resulting dilepton rate is obtained by integrating over the space-time 
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history of the expansion. 

The main result of the previous section is the approximate continuity of the thermal 
emission rate of massive lepton pairs across the phase transition point. Thus the expres- 
sion for the thermal dilepton emission in hadron phase used in hydrodynamical calculations 
should be constructed to ensure this continuity matching. In the narrow mass range, where 



the experimental data on spectral densities is available one can apply Eq. (|T2D with the 
approximate continuity automatically satisfied. To describe the dilepton mass spectrum in 
a wider mass range, as needed to compare with experiment, one needs effective approaches 
like that in |^. Since the binary reactions should be the main source of dileptons at low 
densities and temperatures and since we expect our results to apply best in the vicinity of 
the critical temperature, we will construct an interpolation procedure which at equals the 
quark matter rate and goes over to the binary hadron rate of a low temperatures. 

Below we shall consider the cases of first order and second order deconfinement phase 
transition. The critical temperature in both cases is chosen to be 160 MeV and the freezeout 
temperature 130 MeV. The basic difference between these two cases is the presence or absence 
of a mixed phase. From the arguments of the previous section it is clear, that for the dilepton 
emission the rate should stay equal to the quark-antiquark one for the whole duration of 
mixed phase. 

Let us rewrite a well-known basic formula for the spectrum of dileptons emitted during 



the one-dimensional scaling expansion (see e.g. [Tj], [0) in the form 



dMdy 



.R\ £ drr rX\T{r))^{M,T{r)) , (15) 



where Ra is the radius of the smaller of the two colliding nuclei, r'^q]^'^^ specifies the weight 
for the quark matter and hadron gas thermal rates (lower index) in each of the phases, 
quark matter, mixed or hadron gas (upper index). A head-on collision is assumed. To 
specify completely the evolution of the system one has to know the equations of state for 
both phases and the nature of the phase transition. The entropy of the quark-gluon phase 
will be taken equal to that of the ideal gas of gluons and three lightest quarks 

747r2 1 

^3 _^ 

45 7r2 

A hadron phase will be described by a free resonance gas including all particles with the 
masses up to 2 GeV. The dependence of temperature on proper time T(r) follows from the 



Sq{T) = -^T' + —{UmlTK^imjT) + ?,mlK^{mjT)) . (16) 



Bjorken equation [|T5 



S(T(r))='-^,^-^. (17) 

r uRa dy 
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Let us start with the case when the deconfinement transition from the hadron gas to 
the quark-gluon plasma is of second order. Here the entropy is continuous across the phase 
transition point Tc. We shall assume that at all temperatures the entropy is given by a 



slightly generalized form of an interpolation used in ||17|| ) 

S{T) = (1 - v{T))S,{T) + v{T)S,{T) , (18) 

where 

r^(T) = 1(1 + tanh[(r-r,)/Ar]) (19) 

and AT specifies the width of the transient zone where the entropy density of one phase 
goes over to that of the other phase. We shall take AT = 16 MeV, 10 % of T^. Let us note, 
that the interpolation ([TsD does not imply that the phases are mixed at any temperature, it 
just provides a smooth interpolation between the two known asymptotic expressions. The 
dilepton emission rate is equal to that from quark matter, dB}^^~ / dM"^ , above Tc, i.e. the 
corresponding weights in Eq. ([TsD at T > Tc are = and = 1. Below Tc it is chosen to 
be a combination of the quark matter and hadronic rates: 



dM^ 

The weights in this equation are 



T<Tc 




= r^(T)/(l-r^(T)), 
= l-r,(T) 

ensuring that at T = Tc the overall rate equals that from the quark matter and providing 
an interpolation between the quark rate at Tc and hadron rate at low temperatures. As 
discussed in detailed in the previous section, the necessity of an extra contribution in the 
hadronic phase is due to the necessity of taking into account the effects beyond the binary 
collision approximation, in which the results of were obtained. Otherwise Eq. ([ToD will 
not hold at T = Tc. 

In Fig. we show the dilepton spectra calculated for two different particle densities 
corresponding to head-on Pb-Pb collisions at zero rapidity at SPS {dN/dy = 600, lower 
curves) and RHIC {dN/dy = 1500, upper curves) energies. The dashed curves correspond to 
the binary collision approximation to the hadronic dilepton rate and the solid ones to the rate 
interpolation Eq. (pO|). We see that at SPS energies the contribution from the extra sources 
in hadron phase which are not taken into account in the binary collision approximation 
with a limited number of interacting mesons is visible in the resulting dilepton spectrum 
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Figure 2: Dilepton spectra for SPS (lower curves) and RHIC (upper curves) for second order 
phase transition. The dashed curves correspond to a binary colhsion approximation to the 
hadronic dilepton rate, the solid ones correspond to the rate interpolation Eq. (pO]). 



and brings a multiplicative factor of the order of 1.5. At RHIC energies the signal is already 
almost totally dominated by pairs from quark-antiquark annihilations at the hot early stages 
of the collision (T > Tc). Thus the details of the description of the hadronic electromagnetic 
currents in this mass range, even in the vicinity of Tc, are no longer important. 

Let us now consider the case of a first order phase transition. The basic difference from 
the previous case is that the matter evolves from the quark phase to the hadron one through a 
mixed phase. Since in our approach at T = Tc the rate in hadron phase approximately equals 
that in quark phase, the weights in the mixed phase are = 1 and = 0. This enhances 
the contribution from the mixed phase as compared to the standard scenario in which only 
binary hadron collisions are producing dileptons from a volume occupied by hadronic phase 
throughout the time the mixed phase exists (cf. Fig. 1). In the hadronic phase, T < Tc, the 
rate is taken to be the same superposition of quark and hadron rates, Eqs. (|20|) and (|2T|) , as 
in the case of second order transition. The results of calculation for Pb-Pb collisions at SPS 
and RHIC energies are plotted in Fig. ^. Here the picture is somewhat more complicated 
than in the case of second order transition. The main conclusions can be formulated as 
follows. At SPS energies the system is created at a temperature close to a critical one. 
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Figure 3: Dilepton spectra for SPS (dashed line) and RHIC (solid line) for the first order 
phase transition. 

so the main contribution is coming from the mixed and low-temperature phases. The new 
sources enhance the dilepton rate approximately with a factor of 3 - 4 in the mass range 
under consideration. The contributions originating from the new sources in the mixed and 
the low-temperature phase are of the same order. At RHIC energies the contribution of the 
high temperature phase is already quite significant at M = 2 GeV and dominates at M = 4 
GeV. The new sources enhance the rate with a factor 2 at M ~ 2 Gev and are no longer 
significant at M ~ 4 GeV. 

Let us note, that our main goal in this model calculation is to get a relative enhancement 
of the dilepton rate as compared to the usual calculations exploiting binary collision ap- 
proximation. A well-known drawback of the one-dimensional scaling expansion are the long 
expansion times, ~ 100 fm already at SPS energies. In a more realistic three-dimensional 
treatment the relative time scales for different stages of system's evolution can change, but 
the overall relative enhancement for the contribution from the hadron phase can be expected 
to be more stable. In particular, this should be the case at SPS where the mixed and hadron 
phases are expected to dominate. 
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3 Conclusion 



We have presented a method for taking into account the higher order virial corrections to 
the dilepton emission rate in a hot and dense hadron gas. To go beyond a binary colhsion 
approximation with a hmited number of interacting meson species we argue that the thermal 
emission rates of large mass lepton pairs from the quark matter and dense hadron gas are 
approximately equal in the vicinity of Tc. We then construct an interpolation between 
the hadron rate due to binary collisions at low T and the rate at the critical point T = 
Tc, effectively the quark matter rate. The resulting rates are used in a one- dimensional 
hydrodynamical description of heavy ion collisions at SPS and RHIC energies. The main 
result in the case of a second order deconfinement transition is that while at SPS energies 
the higher order virial corrections to the dilepton rate from the binary hadron collisions 
are still seen in the total thermal dilepton spectrum, at RHIC energies the emission from 
quark matter above dominates in the considered mass range 1.5 GeV < M < 4 GeV. In 
the case of a first order deconfinement phase transition the rate enhancement due to new 
sources is bigger than in the second order case at SPS energies. At RHIC it can still be 
seen at the low-mass end of the mass interval under consideration (M ~ 2 GeV) but is no 
longer important at M = 4 GeV, where the emission from the high temperature QGP phase 
becomes dominant. 

A. Leonidov acknowledges the warm hospitality at the University of Jyvaskyla, where 
a major part of this work was done. His work was partially supported by the Russian 
Foundation for Basic Research under grant 93-02-3815 and by the Academy of Finland 
under grant 27574. 
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